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We analyzed the microlensing of the X-ray and optical emission of the 
lensed quasar PG 1115+080. We find that the effective radius of the X-ray 
\ emission is l-3^ ' 5 dex smaller than that of the optical emission. Viewed as 

a thin disk observed at inclination angle i, the optical accretion disk has a 
scale length, defined by the point where the disk temperature matches the rest 
frame energy of the monitoring band (kT = hc/\ rest with A rest = 0.3/im), of 
\og[{r s opt / cm) a/ cos(i) / 0.5] = 16.6±0.4 . The X-ray emission region (1.4-21.8 keV 
in the rest frame) has an effective half-light radius of log[ri/ 2) x/ cin ] = 15.6^o;g. 
Given an estimated black hole mass of 1.2 x 10 9 M Q , corresponding to a grav- 
itational radius of log[r 9 /cm] = 14.3, the X-ray emission is generated near the 
O ' inner edge of the disk while the optical emission comes from scales slightly larger 
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than those expected for an Eddington-limited thin disk. We find a weak trend 
supporting models with low stellar mass fractions near the lensed images, in 
CN ' mild contradiction to inferences from the stellar velocity dispersion and the time 

00 ■ delays. 
O 
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Introduction 



When iBlaesI (120071 ) recently reviewed the state of accretion disk physics, he found that 
one of the most g l aring probl e ms in even the most sophis ticated accretion disk models 



.e.g. 



Hubeny et al. 



2000 



20011 ; iHirose. Krolik fe Stond 120061 ) is their failure to support a 
hot coron a or to produce X-rays a t all. While there are models for producing the X- 



rays, (e.g. 



Haardt fe Maraschi 



1991 



1993; Maraschi fe Tavecchio 2003; Hirose et al. 2004 



Nayakshin. Cuadra fe Sunyaevll2004l ). they do so on very different physical scales relative to 
the gravitational radius r g = G Mbh/c 2 of the black hol e. For example, the magnetohydro- 
dynamic (MHD) simulations of lHawley fe Balbud (120021 ) predict the dragging of hot ionized 



gas from a jet across the surface of a cooler accretion disk resulting in bremsstrahlung. In 
this model, much of the emission comes from an inner torus with radius r < 20r 5 , but the 
continuum e mission region exten ds to very large radii (r « 200r 9 ). On much smaller scales, 



the model of IHirose et al.l (120041 ) suggests a relativistic MHD accretion disk model in the 



Kerr metric whose inner torus (r < 10r 5 ) supports a large current density capa ble of emitting 



a moderate X-ray flux. The disk-corona model of lHaardt fc Maraschi! ( Il99ll . see also Mer- 
loni 2003) produces X-rays via inverse Compton scattering in a corona which extends over 
much of the o ptical/UV accretion disk, while the "lamp-post" (IMartocchia et al.l 120021 ) and 
"aborted jet" (IGhisellini et al.ll2004l ) models predict a significantly smaller emission structure 



(r <3.0r s ). 

Given their small angular size, few traditional observational constraints can be placed on 
the size of quasar X-ray continuum e mission regions apart from simple a nd often inconclusive 



variability timescale arguments (e.g. IVaughan. Fabian fe Nandrall2003l ). iFabian et al.l (119951 ) 



demonstrated that the broad Fe Ka X-ray emission line in Seyfert 1 spectra is probably 
emitted from the region immediately surrounding the black hole. The width and variability 
of Fe Ka emission has now been m easured in a number of systems to study the innermost 



regions of those accr etion disks (e.g. Ilwasawa et al. 



1999 



Lee et al.l 



999 



Fabian et al 



2002 



Iwasawa et al.ll2004l ). Motivated by this work, lYoung fe Reynolds! ( 120001 ); iBallantyne et al. 
( 120051 ) and others have proposed the use of Constellation- X to measure the size of the Fe 
Ka X-ray reflection region by reverberation mapping. Fortunately, gravitationally lensed 
quasars can be studied on these scales at all wavelengths because the quasar is microlensed 
by the stars in the lens galaxy. The Einstein radius Re of the stars is comparable to 
the expected near-IR sizes of accretion disks, so most disk emission will be significantly 
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microlensed with the amplitude of the variability increasing rapidly for source components 
that are small compared to Re due to the presence of caustic c urves on which the microlensing 
magnification diverges (see the review by IWambsgansd 120061 ) . 



More generally, microlensing variability is a function of the relative tangential velocity v e 
between source, lens and observer, the macroscopic lensing properties of the lens galaxy (the 
convergence k, the stellar surface density fraction and the shear 7) and the relative sizes 
of the source and the source plane projection of the Einstein radius Re of an average mass 
star (M) in the lens galaxy. Since the size of the X-ray emitting region is expected to be much 
smaller than the optical accretion disk, we expect that the ef fects of microlensing will be 



more pronounced at X-ray wavelengths than in the optica l fe.g. iJovanoyic et al.ll2008f). This 



effect has now been observed in many lensed quasars (e.g. iMqrgan et al. 



2002; Dai et al. 2003; Blackburne et al. 20061: Chartas et al. 



was documented specifically in PG 1115+080 by 



2001 



Chartas et al. 



200 



Poolev et all (bood l. 



7UPoolev et al.l l2007h and 



The quadruply-lensed quasar PG 1115+080 was discovered over 25 years ago (jWeymann et al 



19801 ). Since then, it has been the subject of a large number of investigations at multiple 
wavelengths. In particular, the closely separated Al and A2 images bracket a critical line so 
we expect their flux ratio to be approximately unity, but llmpey et al.l (119981 ) and others have 
measured an anomalously low flu x ratio in the optical and NIR (e.g. A2/A1 = 0.64 + 0.02 in 
the if-band). IChiba et all (120051 ) showed that the A2/A1 flux ratio returns to nearly unity 
in the mid-IR (A2/A1 = 0.93 ± 0.06 at 11.7/mi), demonstrat ing that stellar microlen sing is 



the likely cause of the anomaly rather than millilensing (e.g. iKochanek fc Dalall 12004 ) 



Recently, iPooley et al.l (120061 ) conducted a study of the system's anomalo us X-ray flux 



ratios as measured in two Chandra X-Ray Observatory ( Chandra) observations. IPooley et al. 



(120061 ) demonstrated that microlensing is the likely cause of the X-ray flux ratio anomaly 
in PG 1115+080 and qualitatively argued that its X-ray continuum emission region must 
be significantly smaller than its optical accretion disk. In this paper, we combine these 2 
epochs of X-ray data from Chandra with our optical monitoring data to make simultaneous 
measurements of the system's optical and X-ray con t inuum emission regions using the Monte 
Carlo microlensing analysis technique of IKochanek! (120041 ) . 



In § [2j we describe our optical monitoring data and the X-ray flux measurements. In § [31 
we review our microlensing analysis technique and describe its application to PG 1115+080. 
In § |4] we present the results of our calculations and discuss their implications for the sizes 
of the quasar emission regions and the stellar content of the lens galaxy. We assume a flat 
cosmology with Q = 0.3, A = 0.7 and H = 70 km s -1 Mpc -1 . 
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2. Hubble Space Telescope Observations, Chandra and Optical Monitoring Data 



We observed PG 1115+080 in the V- (F555W), I- (F814W) and H- (F160W) bands 
using t he Hubble Space T elescope (HST) for the CfA-Arizona Space Telescope Survey (CAS- 
Falco et al.ll200ll ). The V- and /-band images were taken using th e Wide-Fiel d Plan - 
etary Camera 2 (WFPC2). The if— band images, originally reported in lYoo et al.l (120051 ). 
were taken using the Near-Infrared Camera and Multi-Object Spectrogra ph (NICMOS). W e 
made photometric and astrometric fits to the HST imagery with imfitfits (ILehar et al.ll2000l ). 
using a de Vaucouleurs model for the lens galaxy, an exponential disk model for the quasar 
host gal axy and poin t sourc es for the quasar images. The astrometric fits are consistent with 
those of llmpey et al.l (119981 ). Our HST astrometry and photometry are presented in Table [TJ 



We monitored PG 1115+080 in the i?-band over multi ple seasons with t he SMARTS 

2 the Wisconsin- 



1.3m telescope using the ANDICAM optical/infrared camera (iDePoy et al.ll2003 



I 



the 2.4m 



Yale-Indiana (WIYN) observatory using the WIYN Tip-Tilt Module (WTTM) _ 
telescope at the MD M Observatory using the MDM Eight-K0, Echelle and RET RO C A: 
dMorgan et alihoosh imagers and the 6.5m Magellan Baade telescope using IMACS (jBigelow et al 
19991 ). A d etailed discussion of ou r lensed quasar monitoring data reduction pipeline can 
be found in iKochanek et al.l (120061 ). but we briefly summarize our technique here. We hold 
the lens astrometry fixed to the HST if-band measurements. We treat each quasar image 
as a point source and model the point-spread function with three nested, elliptical Gaussian 
profiles. We measure the flux of each image by comparison to the flux of 5 reference stars 
in the field. We assume that the lens galaxy flux remains constant and fix its value to the 
flux found by minimizing the residuals in a fit to the complete set of measurements from 
each instrument . We suppl emented our opt i cal lig htcurves with ^/-band data published by 
Schechter et al.l (119971 ). The ISchechter et al.l (119971 ) data set does not overlap with any of our 
new monitoring data, so we were unable to correct it for the wavelength difference between 
monitoring bands. Fortunately, this difference is small enough to hav e little effect on the 
resul ts given the expected A 4 / 3 scaling of the optical accretion disk size (IShakura fc Sunvaev 
1973h and our measurement uncertainties. We applied magnitude offsets (e.g. Ofek fc Maoz 
20031 ) to the monitoring data from the other observatories to match the i?-band measure- 



1 lit t p : / /cf a . harvard . edu /castles/ 

2 http:/ /www. astronomy.ohio-state.edu/ANDICAM/ 

3 http: / /www. wiyn.org/wttm/WTTM_manual. html 

4 http:/ /www. astro. Columbia. edu/ arlin/MDM8K/ 

5 http: / /www. astronomy.ohio-state.edu/MDM/RETROCAM 
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merits from SMARTS. 



PG 1115+080 is a particularly challenging system to monitor because the Al and A2 
images are separated by a mere 0'/48. The seeing in our ground-based observations was 
rarely better than l'.'O, so we were forced to sum the flux from ima ges Al and A2. We refer 
to this summed lightcurve as A12 = Al + A2. As documented by iPooley et al.l (120061 ). the 
strongest effect of optical microlensing appears in the A1/A2 flux ratio, so we supplemented 
our lightcurves with 7 epochs of data from the literature in which the Al and A2 images are 



clear ly resolved (ISchechter et al.lll997l ; ICourbin et al .11 199 71 ; llmpey et al.lll998l ; IPooley et al. 
20061 ) . We present our optical monitoring data in Table [31 



We complement our optical ligh tcurves with X-ray f luxes from the two epochs of 0.5 — 
8 keV Chandra imagery publish ed bvlPoolev et al.l (120061 ). although here we used the refined 
flux measurements presented in lPooley et al.l (120071 ). The details of the X-ray data reduction 
and flux ratio calculations are found in those papers. 



3. Microlensing Models 



Microlen sing; statistics are strongly influe nced by the presence of smoothly distributed 
dark matter (ISchechter &: Wambsgansd 120021 ) . typically parameterized as the ratio of the 
stellar surface density to the total surface density We considered a range of possible 

stella r mass fractions in our calculations. We used the GRAVLENS software package (IKeeton 
20011 ) to generate a series of ten mode ls that match the HST astrometry and reproduce the 
mid-infrared (11.7/im) flux ratios fromlchiba et al. ( 2005 ). Each model consists of concentric 
de Vaucouleurs and NFW (jNavarro. Frenk fc White! Il996l ) profiles, and we vary the mass in 
the de Vaucouleurs component over the range 0.1 < fu/L < 1-0 in steps of A/m/l = 0.1, 
where Jm/l — 1-0 represents a constant mass-to-light ratio (de Vaucouleurs) model with 
no da rk matter halo. T a ble [2] summarizes the microlensing parameters as a function of 
/m/l- iTreu fc Koopmansl (120021 ) found that th e best fit to the system's large stellar velocity 
dispersion (cr* = 281 ±25 km s" 1 . lTonrylll998l ) is provided by a steep mass profile p cx r -2 35 , 
implying a large stella r mass component, and for Ho = 72 km s _1 Mpc -1 , the best fit to the 



Schechter et al.l (119971 ) time delays is provided by the / 



M/L 



0.8 model. 



We generated a set of microlensing magnification patterns at each image l ocation for each 
of the 10 macroscopic mass models using a variant of the ray-shooting method (jSchneider et al 



19921 . see iKochanekl 120041 for the details of our technique). The patterns are 8192 x 8192 



images of the source-plane projection of the magnification patterns from an ensemble of 
typical lens galaxy stars at each image location. We approximated the Galactic stellar mass 
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function of iGouldl (120001 ) as a power law, assuming dN(M)/dM oc M~ 13 with a dynamic 
range in mass of a factor of 50. The mean stellar mass in the lens galaxy (M) is initially un- 
known, so magnification patterns are produced in units of the Einstein radius with an outer 
scale of 20 R E . For PG 1115+080, the Einstein radius is R E = 6.6 x 1O 16 (M/M ) 1 / 2 C m. 
To convert to physical units, all results are eventually scaled by some factor of (M/M Q ). 
Henceforth, quantities in Einstein units will be given the "hat" accent to distinguish them 
from quantities in physical units. So the physical source size r s is related to the scaled source 
size r s by r s = f s (M/M ) 1//2 , and the physical effective velocity v e is related to the scaled 
velocity v e by v e = v e (M/M Q ) 1/2 . 

In order to eliminate the quasar's intrinsic vari ability, we shifted the optical light curves 
by the measured time delays (jSchechter et al.lll997l ). so that any remaining variability in the 
flux ratios must be attributed to microlensing. It is impossible to offset the sparse X-ray 
flux measurements by the time delays, so we assume that X-ray flux ratios can be treated 
as simultaneous in a statistical sense. The time delay between the Al and A2 images is 
less than one day, so there was no need to apply a time delay correction to the 7 epochs of 
individually resolved Al and A2 data. 

As described in detail by lKochanek! (12004 ) (see also lMorgan et al. 1 120071 : iPoindexter. Morgan fc Kochan 



20081 ). our Monte Carlo microlensing analysis searches for trajectories across the magnifi- 
cation pattern that fit the observed li ght curves. We used a th in accretion disk surface 
brightness profile for the source model (jShakura fc Sunyaevlll973l ) with 



I(R) oc <^ exp (R/r, 



\3/4 



1 



where the scale radius r s is the radius at which the disk temperature matches the rest-fame 
wavelength of our monitoring band, kT = hc(l + z s )/X f >s - We neglect the central hole in 
the emission profile, the effect of which is negligible at optical wavelengths. Microlensing 
primarily depends on the projected area of the source while the true scale lengths also depend 
on the shape of the source and its inclination. We will refer to a radius where we have ignored 
the shape and inclination of the source as an "effective" radius that defines a projected 
For a thin disk, the effective radius is related to the source scale length by 
r 2 s cos i where i is the inclination angle. The X-ray emissio n presumably has a different 



area irr 



'eff 

emission profile and shape. Fortunately. iMortonson et al.l (J2005|) demonstrated that the half- 
light radius measured with microlensing is essentially independent of the surface brightness 
profile, so we will characterize the X-ray emission by the effective half-light radius. For our 
thin disk model, the half-light radius is related to the disk scale length by R\ii = 2A4r s . 
In summary, to compare the sizes of the optical and X-ray emitting regions we will use the 
ratio of the effective radii r opt /rx, to characterize the optical emission we will use the thin 
disk scale length r St0pt and an inclination angle cos i, and for the X-ray emission we will use 
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the effective (i.e. no shape corrections) half-light radius T\iix- 



We generated 8 trial magnification patterns for each of the 10 macroscopic mass models. 
For each trial and model we produced 50,000 trial light curves for a 16 x 21 grid of X-ray and 
optical source sizes. These source sizes r s are scaled sizes that depend on the microlens mass 
r s = f s (M/M ) 1 / 2 . We used logarithmic grids spanning the region producing acceptable 
fits with a grid spacing of 0.2 dex. In total there were 4 x 10 6 trial light curves for each 
combination of X-ray and optical source sizes. When assessing the quality of our fits to 
the observed flux ratios, we allowed for only 0.1 mag of syst ematic unce r tainty in the flux 
ratios of the macro models because the mid-IR flux ratios of Chiba et al. J2005h are a close 
approximation to the intrinsic flux ratios in this system. In selecting trial light curves we 
gave equal statistical weight to the optical data where Al and A2 could not be separately 
measured, optical data where Al and A2 could be separately measured, and the X-ray data 
so that we would isolate trials with reasonable fits to all three classes of data. The final 
goodness of fit was evaluated with a x 2 At to the light curves where all data have their true 
statistical weights and we discard all fits with \ 2 /Nd„f > 4.0 as they make no significant 
contribution to the final Bayesian integrals (see iKochanek! 120041 ) . Figures [H and [2] show two 
examples of good fits to the data. The stronger flux anomalies in the X-ray data force the 
X-ray source to be more compact than the optical, leading to the much larger variability 
predicted for the X-ray bands relative to the optical. 

To convert the results to physical units, we assume a prior on the mean mass of the 
microlenses (M) and the transverse velocity between source, lens and observer v e . For the 
mean stellar mass prior, we assume O.1M < (M) < l.OM . We model the effective velocity 
of the system with three components. We set t he velocity of the o bserver v Q = 94 km s _1 to 



be the projection of the CMB dipole velocity (IKogut et al.lll9931 ) onto the lens plane. We 
calculate a one-dimensional stellar velocity dispersion in the lens galaxy of cr* = 220 km s _1 
based on the Einstein radius of its macroscopic mass model and we assume a lens ga laxy 
peculiar velocity dispersion of a v = 235/(1 + z\) km s _1 = 179 km s _1 (jKochanek! 120041 ) . In 
general these two priors give similar physical size estimates if applied separately because the 
sizes depend only weakly on the mass scale (M). Microlensing depends only on the size and 
velocity of the source in Einstein units, r s = r s /(M/M ) 1//2 and v e = v e / (M / Mq) 1 ' 2 , and 
a given level of variability can be produced either by moving a small source slowly (both 
r s and v e small) or a large source rapidly (both f s and v e large) with (roughly) r s oc v e . 
But the mass scale implied by a given Einstein velocity scales as (M/M Q ) = (v e /v e ) 2 so the 
dependence of the physical scale on the mass essentially cancels given som e knowledge of the 



physical velocity v e , with r s = r s {M/M & ) 1 / 2 oc (M)° fsee lKochanekH2004f ). For PG 1115+080 
the poor temporal overlap of the optical and X-ray light curves means that the differences 
between using the priors separately are larger than we have found for most other lenses (e.g. 
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Morgan et al.ll2007l ; iPoindexter. Morgan &: Kochanekll2008l ). so we restricted our analysis to 
using the priors jointly. 



Results and Discussion 



Figure [3] shows the ratio of the effective radii of the optical and X-ray sources where 
the effective radius should be viewed as the square root of the projected source area. The 
advantage of the size ratio is that it has no direct dependence on the mass of the microlenses 
(in the sense that f opt /fx = r opt /rx)- We find that \og[r opt /rx] = 1.3~t Q 'l. Figure H] shows 
the estimates for the physical sizes, where we show the inclination corrected disk scale length 
for the optical source and the effective half-light radius for the X-ray source. Recall from 
§ [3] that the disk scale length is the point where the temperature equals the photon energy 
kT = hc/X res t and that the effective half-light radius has no shape or inclination corrections. 
Thus, at 0.3/mi or T = 4.8 x 10 4 K, the disk scale length is log[(r Si0p t/cm) -^/cos(z)/0.5] = 
16.6 ± 0.4, and in the (rest-frame) 1.4-21.8 keV band the effective X-ray half light radius is 
log(r 1/2)X /cm) = 15.6lg;jj. 

We can compare these size esti mates to th e oretic al expectations given the estimated 
black hole mass of 1.2 x 10 9 M Q from lPeng et all (120061 ) based on the quasar luminosity and 
the width of the Mg II (A2798A) emission line, where b lack hole mass estimates using this 



techn ique have a typical uncertainties of ~ 0.3 dex (see iMcLure fc Jarvid 120021 ; iPeng et al. 



20061 ). Fig. H] shows the gravitational radius r g = 1.8 x 10 14 cm for this mass, which is the 



innermost stable circular orbit for a maximally rotating Kerr black hole. For reference, we 
also plot the innermost stable circular orbit for a Schwarzschild black hole at 6r g . The optical 
emission comes from well outside the inner edge of the di sk (~ 220r p ), justifying ou r neglect 
of the inner edge of the disk in Eqn. [TJ Thin disk theory (IShakura &: Sunyaevlll973l ) predicts 
that the optical size for a face-on quasar radiating with 10% efficiency at the Eddington 
limit should be log [ r.,^/ cm] = 15.6, and our black hole mass/accretion disk size scaling 
from iMorgan et al.l (120071 ) predicts a scale radius of lo g [(r,,^ /cm) ^/ c os(i) / 0.5] = 15.8 ±0.3. 
Our current result is marginally inconsistent with the IMorgan et al.l (120071 ) black hole mass 
- accretion disk size scaling, but the small discrepancy could be explained by adjusting the 
unknown dis k inclination angle. Our result is also somewhat larger than the theoretical thin 
disk size (see Morgan et al.ll2007l . for a detailed comparison of microlensing disk size estimates 
in 11 systems). The far more compact X-ray emission comes from a region very close to the 
inner disk edge, with an effective half-light radius of ~ 20r„. Th is result seems to favor 



models with a small o r moderately-sized em i ssion structures (e .g. lHawley &: Balbusl 12002 



Martocchia et al.l 120021 ; iGhisellini et al.ll2004l ; iHirose et al.l 120041 ) and disfavor the standard 
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disk-corona model ( Haardt fc Maraschi 1991 ; Merloni 2003 ) 



We also obtain some information on the structure of the lens galaxy, as illustrated 
in Figure [5] where we show our estimates of the stellar fraction J'm/l- We do not find 
a strong peak in the fu/L distribution, but we do detect a weak trend favoring models 
with lower f M /r , and a low stellar surf ace density (Fig . EJ), as w e woul d expe ct. Both the 
Schechter et al.l (119971) time delays (s ee lKochanekll2002l ) and the iTonryl (119981 ) velocity dis- 



Treu fc Koopmand |2002| ) requ ire mass models with little dark matter near the 



persion (see 

radius of the lensed images (although see iRomanowsky fc Kochanekl Il999l for examples of 
dynamical models consistent with both a significant dark matter halo and the high velocity 
dispersion). While Figure [5] is not conclusive, it is probable that with better X-ray light 
curves we will be able to measure Jmil and either confirm or reject the time delay and 
velocity dispersion measurements. 

We are expanding our analyses to include all 10 lensed quasars with archival X-ray 



data (see lPoolev et al.l 120071 ) as well as three systems (HE 1104-1805, RX J1131-1231 and 



Q 2237+0305) where we have obtained X-ray light curves. For many of these systems 
(RX J1131 1231, Q 2237+0305, WFI J2033-4723, SDSS 09 24+0219 and H 141 3 +117) we 



have optical light curves comparable to those used here (see iMorgan et al.l 120061 ; iDai et al. 



20081 ). but for the remainder we will have to rely on sparse, archival optical data. The 
main challenge we face is that the computational intensity of modeling the two bands si- 
multaneously is a significant bottleneck for completing the analyses. Nonetheless, we see 
no fundamental barriers to complem enting our correlations between optical disk size and 
black hole mass (IMorgan et al.l 120071 ) with their X-ray equivalents. The pattern suggested 
by PG 1115+080 is that the X-ray continuum emission region tracks the inner edge of the 
accretion disk. We hope to determine if this result is universal. 
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Table I. HST Astrometry and Photometry of PG11 15+080 



Component 


Astrometry 




Photometry 






ARA 


ADcc 


H=F160W 


I=F814W 


V=F555W 


Al 


l'/328 ± 0C003 


-2'.'034 ± (X'003 


15.71 ± 0.02 


16.42 ±0.02 


16.90 ±0.11 


A2 


1'.'477 ± 0'.'004 


-1'.'576 ± 0'.'003 


16.21 ± 0.02 


16.85 ±0.01 


17.62 ±0.09 


B 


-0'.'341 ± 0'.'003 


-1'.'961 ± (X'003 


17.70 ± 0.02 


18.37 ±0.01 


18.95 ±0.12 


C 


= 


= 


17.23 ± 0.03 


17.91 ±0.02 


18.39 ±0.06 


G 


0"381 ± 0"003 


-l'/344 ± 0'/003 


16.66 ± 0.04 


18.92 ±0.02 


20.74 ±0.03 



Table 2. PG1 115+080 Lens Galaxy Mass Models 



Im/l 




Convergence k 






Shear 7 






K 








Al 


A2 


B 


C 


Al 


A2 


B 


C 


Al 


A2 


B 


C 


0.1 


0.75 


0.76 


0.84 


0.69 


0.22 


0.27 


0.32 


0.18 


0.021 


0.023 


0.031 


0.016 


0.2 


0.69 


0.71 


0.79 


0.64 


0.27 


0.33 


0.41 


0.20 


0.046 


0.050 


0.066 


0.034 


0.3 


0.64 


0.66 


0.74 


0.59 


0.31 


0.39 


0.49 


0.23 


0.075 


0.081 


0.105 


0.055 


0.4 


0.59 


0.61 


0.69 


0.54 


0.36 


0.44 


0.58 


0.25 


0.11 


0.12 


0.15 


0.08 


0.5 


0.54 


0.56 


0.63 


0.49 


0.40 


0.50 


0.66 


0.28 


0.15 


0.16 


0.20 


0.11 


0.6 


0.49 


0.51 


0.59 


0.44 


0.44 


0.56 


0.74 


0.30 


0.20 


0.21 


0.25 


0.14 


0.7 


0.45 


0.46 


0.54 


0.39 


0.48 


0.61 


0.81 


0.32 


0.25 


0.27 


0.31 


0.18 


0.8 


0.40 


0.41 


0.49 


0.34 


0.52 


0.66 


0.89 


0.35 


0.32 


0.34 


0.39 


0.24 


0.9 


0.36 


0.38 


0.44 


0.30 


0.55 


0.71 


0.96 


0.36 


0.40 


0.42 


0.46 


0.29 


1.0 


0.31 


0.33 


0.40 


0.26 


0.60 


0.76 


1.03 


0.39 


0.51 


0.54 


0.57 


0.39 



Note. — Convergence k, shear 7 and the fraction of the total surface density composed of stars re* /ft 
at each image location for the series of macroscopic mass models where /m/l = 1-0 corresponds to a 
constant mass-to-light ratio model for the lens galaxy. 
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Fig. 1. — Examples of good fits to the observed flux ratios. Top panel: We fit the 7 epochs 
of data with a resolved A1/A2 flux ratio individually. The error bars are the data, and the 
black curve is the fit. Middle panel: Dotted red and dashed blue curves are best fits to the 
A12/B and AY1/C flux ratios, respectively. Data are plotted with error bars in the same 
color scheme. The AY2/B and AYljC flux ratios varied little over the last decade. Bottom 
panel: The observed A1/A2, Al/B and Al/C X-ray flux ratios are plotted using solid black 
triangles, solid red circles and open blue triangles, respectively. The best fits to the observed 
A1/A2, Al/B, and Al/C flux ratios are plotted using solid black, dotted red and dashed 
blue curves, respectively. 
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Fig. 2. — A second solution plotted as in Fig. [TJ 
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X — ray size 
( r i/2,x) 




Optical size 
(r t ) 

v s.opt/ 



source size 



Fig. 4. — Probability distributions for the effective X-ray half light radius ri/2,x (top) and 
optical thin disk scale radius r s opt (bottom). For the thin disk, we assumed an inclination 
of co si = 1/2. Given the black hole mass estimate of 1.2 x 10 9 M for PG 1115+080 



from iPeng et all (120061 ). the solid vertical lines indicate the innermost stable circular orbit 



r g = GMbh/c 2 for a maximally rotating Kerr black hole and the innermost stable circular 
orbit for a Schwarzschild black hole at 6r g . The dashed vertical line is the prediction of thin 
disk theory for the scale radius at 0.26/zm for an Eddington-limited accretion disk radiating 
at 10% efficiency. 
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0.2 0.4 0.6 0.8 1 

macroscopic mass model f M/L 

Fig. 5. — Probability distribution for Jm/l, the fraction of the lens galaxy mass in the 
constant M/L ratio (de Vaucouleurs) component. Jmil can be related to the stellar surface 
density fraction k^/k at each image location using the data in Table El The Jmil value 
implied by the time delays is plotted with a solid vertical line. 
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Table 3. PG1115+080 Optical Light Curves 



HJD 


X 2 /N dof 


T v . . ....... A 1 I AO 

Images J\L-\-l\z 




Image B 




Image C 


^otarsy 


Source 


2994 


.927 


1 


.00 


-1.095 ± 0.006 


1 


.192 ± 0.010 





.896 ± 0.008 


0.019 ± 0.004 


MDM 


3038 


.714 


1 


23 


-1.071 ± 0.008 


1 


.319 ± 0.012 





971 ± 0.010 


-0.014 ± 0.005 


SMARTS 


3045 


.790 


1 


.18 


-1.078 ± 0.014 


1 


.323 ± 0.018 


1 


.064 ± 0.024 


-0.046 ± 0.005 


SMARTS 


3047 


.003 


1 


.69 


-1.078 ± 0.004 


1 


.336 ± 0.007 





.938 ± 0.007 


0.148 ± 0.004 


WIYN 


3055 


.690 


1 


46 


-1.077 ± 0.009 


1 


.321 ±0.012 





973 ± 0.010 


-0.003 ± 0.004 


SMARTS 


3055 


.753 


1 


.09 


-1.080 ± 0.004 


1 


.314 ± 0.008 





.950 ± 0.007 


0.068 ± 0.004 


MAGELLAN 


3063 


.677 


1 


.47 


-1.082 ± 0.008 


1 


.362 ±0.011 





.968 ± 0.009 


0.004 ± 0.004 


SMARTS 


3064 


.824 


1 


55 


-1.077 ± 0.005 


1 


.235 ±0.008 


0. 


900 ± 0.007 


0.066 ± 0.004 


MDM 


3080 


.640 


1. 


.90 


-1.083 ± 0.008 


1 


.374 ±0.011 


0. 


.993 ± 0.009 


0.005 ± 0.004 


SMARTS 


3091 


.573 


1 


.14 


-1.078 ± 0.009 


1 


.309 ±0.012 





938 ± 0.010 


-0.006 ± 0.004 


SMARTS 


3101 


.647 


1 


11 


-1.102 ± 0.011 


1 


.377 ± 0.017 





.918 ± 0.012 


-0.041 ± 0.005 


SMARTS 


3104 


.646 


2 


.85 


-1.105 ± 0.004 


1 


.297 ± 0.007 





.917 ± 0.007 


0.172 ± 0.004 


WIYN 


3108 


.540 


3 


46 


-1.119 ± 0.008 


1 


.289 ±0.012 


0. 


964 ±0.010 


0.000 ± 0.004 


SMARTS 


3116 


.582 


1. 


02 


-1.119 ± 0.009 


1 


.337 ±0.013 





953 ± 0.010 


-0.002 ± 0.004 


SMARTS 


3132 


.534 


1 


.07 


-1.107 ± 0.009 


1 


.300 ± 0.012 





961 ± 0.010 


-0.021 ± 0.005 


SMARTS 


3136 


.717 


5 


.84 


-1.132 ± 0.005 


1 


.255 ± 0.008 





893 ± 0.007 


0.062 ± 0.004 


MDM 


3138 


.476 


1 


34 


-1.093 ± 0.009 


1 


.322 ± 0.013 


1 


.004 ± 0.011 


-0.016 ± 0.005 


SMARTS 


3359 


787 


2 


.12 


-0.929 ± 0.011 


1 


.350 ± 0.014 


1 


.085 ± 0.011 


-0.006 ± 0.004 


SMARTS 


3368 


759 





.66 


-0.888 ± 0.023 


1 


.445 ± 0.034 


1 


.251 ± 0.026 


-0.021 ± 0.005 


SMARTS 


3393 


.988 


3 


.54 


-0.867 ± 0.009 


1 


.690 ± 0.012 


1 


245 ± 0.009 


0.032 ± 0.004 


MDM 


3394 


.754 





.80 


-0.870 ± 0.015 


1 


.531 ±0.019 


1 


.117 ± 0.031 


-0.040 ± 0.005 


SMARTS 


3403 


.785 


1. 


.00 


-0.850 ± 0.011 


1 


.555 ± 0.014 


1 


.199 ± 0.011 


-0.011 ± 0.004 


SMARTS 


3413 


792 


2 


84 


-0.850 ± 0.009 


1 


.694 ±0.014 


1 


286 ± 0.011 


0.001 ± 0.004 


SMARTS 


3417 


.760 


1 


.55 


-0.832 ± 0.009 


1 


.637 ± 0.013 


1 


307 ± 0.015 


-0.019 ± 0.004 


SMARTS 


3424 


727 





.60 


-0.818 ± 0.022 


1 


.672 ± 0.035 


1 


.216 ± 0.038 


-0.026 ± 0.005 


SMARTS 


3428 


.737 


1 


19 


-0.791 ± 0.014 


1 


.709 ± 0.024 


1 


.354 ± 0.031 


-0.042 ± 0.005 


SMARTS 


3431 


786 


1 


51 


-0.836 ± 0.014 


1 


.772 ± 0.020 


1 


282 ± 0.014 


-0.030 ± 0.005 


SMARTS 


3433 


.726 


1 


.44 


-0.832 ± 0.010 


1 


.689 ±0.015 


1 


.230 ± 0.011 


-0.009 ± 0.004 


SMARTS 


3435 


.754 


1 


.71 


-0.830 ± 0.008 


1 


.634 ± 0.013 


1 


.224 ± 0.010 


0.006 ± 0.004 


SMARTS 


3442 


.691 


2. 


05 


—0.833 ± 0.011 


1 


.607 ± 0.015 


1 


245 ± 0.011 


—0.010 ± 0.004 


SMARTS 


3445 


.734 


2 


.31 


-0.845 ± 0.009 


1 


.696 ±0.015 


1 


.254 ±0.011 


-0.005 ± 0.004 


SMARTS 


3447 


.633 


0. 


.82 


-0.833 ± 0.013 


1 


.630 ±0.017 


1 


.199 ± 0.012 


-0.018 ± 0.005 


SMARTS 


3449 


.741 


1 


39 


-0.823 ± 0.010 


1 


.630 ±0.014 


1 


232 ± 0.011 


-0.013 ± 0.004 


SMARTS 


3458 


.682 


0. 


.96 


-0.819 ± 0.012 


1 


.675 ±0.017 


1 


.264 ±0.013 


-0.032 ± 0.005 


SMARTS 


3459 


.599 





.81 


-0.818 ± 0.015 


1 


.604 ± 0.020 


1 


.244 ± 0.014 


-0.030 ± 0.005 


SMARTS 


3461 


.676 


0. 


86 


-0.825 ± 0.012 


1 


.641 ±0.016 


1 


211 ± 0.012 


-0.010 ±0.005 


SMARTS 


3468 


.600 


2. 


00 


-0.828 ± 0.009 


1 


.630 ±0.014 


1 


253 ± 0.011 


-0.001 ± 0.004 


SMARTS 


3470 


.608 


2 


.68 


-0.837 ±0.011 


1 


.543 ±0.015 


1 


.324 ±0.012 


-0.006 ± 0.004 


SMARTS 


3478 


.556 


0. 


83 


-0.810 ±0.014 


1 


.595 ±0.018 


1 


214 ±0.013 


-0.033 ± 0.005 


SMARTS 


3483 


.539 


1 


05 


-0.793 ± 0.017 


1 


.603 ± 0.026 


1 


677 ± 0.030 


-0.037 ± 0.005 


SMARTS 


3490 


.600 


1 


.49 


-0.822 ± 0.011 


1 


.653 ±0.015 


1 


.256 ± 0.011 


-0.013 ± 0.004 


SMARTS 


3491 


.594 


1 


.34 


-0.808 ± 0.011 


1 


.688 ±0.015 


1 


263 ± 0.011 


-0.011 ± 0.004 


SMARTS 


3523 


.497 


0. 


97 


-0.754 ±0.012 


1 


.703 ±0.017 


1 


270 ±0.012 


-0.017 ±0.005 


SMARTS 


3527 


.474 


1 


.20 


-0.761 ± 0.011 


1 


.753 ±0.016 


1 


.318 ± 0.012 


-0.014 ±0.005 


SMARTS 


3718 


.838 


0. 


.92 


-0.787 ±0.015 


1 


.562 ±0.018 


1 


.174 ±0.013 


-0.024 ±0.005 


SMARTS 


3737 


.018 


4. 


.71 


-0.885 ± 0.004 


1 


.434 ± 0.007 


1 


127 ±0.007 


0.099 ± 0.004 


MDM 



Table 3 — Continued 



HJD X 2 l N dof Images A1+A2 Image B Image C (Stars) Source 



3750 


.790 


0.68 


-0.795 ± 0. 


022 


1.663 ±0. 


035 


1.331 ± 0. 


026 


-0.030 ± 0. 


005 


SMARTS 


3755 


.748 


0.63 


-0.784 ± 0. 


023 


1.613 ±0. 


037 


1.186 ± 0. 


024 


-0.020 ±0. 


005 


SMARTS 


3761 


.765 


1.09 


-0.820 ±0. 


010 


1.579 ±0. 


013 


1.183 ± 0. 


010 


0.005 ± 0. 


004 


SMARTS 


3772 


.763 


1.71 


-0.828 ± 0. 


.009 


1.547 ±0 


.013 


1.175 ± 0. 


010 


0.007 ±0. 


004 


SMARTS 


3792 


.766 


2.47 


-0.845 ± 0. 


Oil 


1.706 ±0. 


.017 


1.240 ±0. 


012 


-0.006 ± 0. 


004 


SMARTS 


3815 


.678 


1.48 


-0.851 ±0 


.010 


1.721 ±0. 


.017 


1.116 ± 0. 


.012 


-0.012 ±0. 


005 


SMARTS 


3823 


.631 


2.25 


-0.885 ± 0. 


.013 


1.646 ±0 


018 


1.228 ± 0. 


.013 


-0.016 ± 0. 


.005 


SMARTS 


3855 


.526 


0.59 


-0.883 ± 0. 


030 


1.522 ±0. 


.058 


0.996 ± 


051 


0.023 ± 0. 


.006 


SMARTS 


3884 


.450 


1.00 


-0.908 ± 0. 


.019 


1.502 ±0. 


.024 


1.118 ± 


.024 


-0.033 ± 


.005 


SMARTS 


3893 


.475 


0.64 


-0.921 ± 0. 


.025 


1.655 ± 


.035 


1.130 ±0 


.022 


-0.032 ± 


.005 


SMARTS 


3900 


.490 


0.94 


-0.919 ± 0. 


.012 


1.548 ±0 


017 


1.159 ± 


.012 


-0.023 ± 0. 


.005 


SMARTS 



Note. — HJD is the Heliocentric Julian Day -2450000 days. The goodness of fit of the image, \ 2 /^dofj 
is used to rescale the formal uncertainties by a factor of (x 2 l^dof) 1 ^ 2 ■ The Image columns give the 
magnitudes of the quasar images relative to the comparison stars, and the Images A1+A2 column is the 
sum of the flux from both images. The (Stars) column gives the mean magnitude of the standard stars for 
that epoch relative to their mean for all epochs. 



